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Abstract 

Although  there  have  been  no  major  revisions  of  the  National  Bureau  of  Standards  Handbook  of  Mathe- 
matical Functions  [1]  since  its  first-  publication  in  1964,  it  continues  to  be  widely  used  by  the  mathematical 
and  scientific  community.  As  a result,  NIST  is  engaged  in  a large  scale  project  to  update  and  expand 
the  handbook  and  disseminate  it  in  digital  format  on  the  World  Wide  Web.  This  electronic  handbook 
is  being  called  the  Digital  Library  of  Mathematical  Functions  (DLMF).  A key  feature  of  the  DLMF 
will  be  dynamic  3D  visualizations  of  special  functions.  This  paper  discusses  some  of  the  visualization 
results  obtained  to  date  using  VRML,  the  Virtual  Reality  Modeling  Language,  a standard  file  format  for 
describing  the  behavior  and  geometry  of  a 3D  virtual  world. 

Keywords:  VRML,  Virtual  Reality  Modeling  Language,  visualization,  3D  graphics,  clipping,  special 
functions 


Introduction 

One  of  the  most  well  known  publications  to  come  out  of  the  National  Bureau  of  Standards,  the  prede- 
cessor organization  of  the  National  Institute  of  Standards  and  Technology  (NIST),  is  the  Handbook  of 
Mathematical  Functions  [1],  published  in  1964.  LInchanged  since  the  tenth  printing  in  1972,  it  continues 
to  be  widely  sold  by  the  US  Government  Printing  Office,  Dover,  and  many  other  commercial  publishers. 
The  Science  Citation  Index  indicates  that  the  handbook  was  cited  over  26,000  times  during  the  period 
1974-1997  in  more  than  1000  journals,  and  the  number  of  citations  continues  to  increase  with  time. 

The  continued  interest  in  the  handbook  plus  such  factors  as 

• the  clear  advantages  of  electronic  media  for  the  construction  and  communication  of  ideas  in  technical 
fields; 

• advances  in  basic  mathematical  and  computational  techniques  associated  with  the  classical  special 
functions  of  the  mathematical  and  physical  sciences;  and 

• the  identification  of  new  functions  having  widespread  importance  in  emerging  applications 
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have  led  NIST  to  embark  on  a massive  project  to  update  and  expand  the  current  handbook  and  dis- 
seminate it  in  digital  format  on  the  World  Wide  Web;  see  Lozier  [2]  for  an  early  description  and 
http://math.nist.gov/DigitalMathLib  for  current  information  on  this  new  project.  The  new  entity,  which 
is  being  called  the  Digital  Library  of  Mathematical  Functions  (DLMF),  will  make  full  irse  of  advanced 
communications  and  computational  resources.  A key  feature  of  the  DLMF  will  be  dynamic  3D  visualiza- 
tions of  special  functions  that  allow  a user  to  conduct  interactive  explorations  of  the  relationship  between 
a function’s  mathematical  or  numerical  properties  and  its  graphical  representations.  The  Virtual  Reality 
Modeling  Language  (VRML)  is  a standard  3D  file  format  for  describing  the  behavior  and  geometry  of 
a 3D  virtual  world,  or  scene.  Its  accessibility  on  the  Internet  and  interactive  capabilities  make  it  an 
ideal  candidate  for  this  development  work.  This  paper  describes  techniques  we  are  using  to  generate 
VRML  hies  that  not  only  display  3D  graphs,  but  also  allow  the  user  to  gain  a deeper  understanding  of 
the  function  by  manipulating  cutting  planes  perpendicular  to  the  major  axes  in  order  to  examine  the  2D 
projection  of  the  function  at  various  locations. 


3D  visualization  of  special  functions 

Like  the  original  handbook,  the  DLMF  is  designed  primarily  for  the  use  of  scientists.  A secondary, 
but  important  goal  is  to  reach  a much  broader  audience  by  making  aspects  of  the  DLMF  accessible  to 
educators  and  students.  An  obvious  way  to  support  these  dual  goals  is  to  create  3D  visualizations  that 
are  both  exciting  and  informative.  Fortunately,  the  graphical  representations  of  many  special  functions 
are  so  complex  and  interesting  that  by  designing  visualizations  that  illustrate  the  features  of  interest  to 
scientists  we  automatically  produce  displays  that  grab  the  attention  of  less  technically  oriented  viewers. 
In  any  case,  the  development  of  the  display  requires  close  coordination  with  an  expert  in  the  field  of 
special  functions.  Currently,  we  are  concentrating  on  visualizations  for  the  chapter  on  Airy  functions, 
written  by  Prof.  Frank  Olver,  one  of  the  authors  of  the  original  handbook.  The  Airy  functions,  Ai  and 
Bi,  occur  in  quantum  mechanics,  in  the  study  of  wave  diffraction,  electromagnetism,  and  other  areas  of 
physics  and  engineering,  and  arise  as  solutions  of  the  second  order  differential  equation 

d2w 

~d^  = ZW- 

To  obtain  reliable  data  for  the  visualizations  we  used  a double  precision  Fortran  routine  for  the 
calculation  of  Airy  functions  written  by  D.E.Amos  [3].  We  wrote  a C program  that  accepted  this  input 
data  and  generated  a VRML  file  as  output.  In  many  cases  we  found  that  the  graphs  had  to  be  scaled 
very  carefully  in  order  to  make  interesting  features  visible.  However,  for  some  functions,  the  variation  in 
values  over  the  domain  was  so  extreme  that  simply  adjusting  the  scaling  was  not  enough.  A 3D  clipping 
algorithm  could  potentially  help  with  the  scaling,  but  to  date  we  have  been  unsuccessful  in  finding  suitable 
routines  in  the  literature  or  in  available  packages.  We  first  tried  resetting  values  above  a certain  height  to 
the  same  constant  as  was  done  by  Thompson  [4],  but  that  produced  the  misleading  “table”  effect  shown 
in  Figure  1.  We  also  tried  suppressing  the  plotting  of  points  where  the  function  value  was  greater  than 
a specified  constant,  but  that  produced  plots  with  jagged  edges  that  were  equally  misleading.  Finally, 
we  decided  to  use  information  from  a contour  plot  of  the  function  to  restrict  the  domain  to  points  where 
the  function  values  were  less  than  or  equal  to  a specified  constant.  The  contour  information  was  used  to 
construct  the  boundary  of  the  domain,  and  a boundary  fitted  mesh  was  then  placed  over  the  domain  as 
shown  in  Figure  2.  By  computing  the  Airy  function  only  at  values  on  the  mesh,  we  obtained  a smoothly 
clipped  surface  plot  as  shown  in  Figure  3. 
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Figure  1:  |Bi'(,z)|,  Modulus  of  the  derivative  of  Bi(2). 


Figure  2:  Contour  mesh. 
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Figure  3:  Clipped  version  of  |Bi'(z)|. 


Intersection  of  3D  surfaces  with  cutting  planes 

While  a VRML  browser  gives  the  user  the  flexibility  to  examine  a 3D  display  from  an  arbitrary  direction, 
we  also  wnnted  to  give  the  user  a greater  ability  to  interact  with  the  visualization  of  the  functions.  DLMF 
project  members  noted  that  one  very  useful  feature  would  be  the  ability  to  study  the  changes  that  occur 
in  the  intersection  of  a plane  with  a surface  as  the  plane  is  moved  through  the  surface.  Using  VRML  we 
were  able  to  create  files  that  would  generate  cutting  planes  through  a 3D  surface.  However,  currently, 
the  cutting  planes  are  limited  to  planes  perpendicular  to  the  X and  Y coordinate  axes. 

We  defined  a VRML  PROTO  called  Cutplane  to  display  the  plane  and  compute  the  intersection  of 
the  plane  with  the  surface  data.  In  VRML  a PROTO  is  a user  defined  reusable  component.  Cutplane 
searches  the  data  to  determine  which  surface  points  are  closest  to  the  specified  plane  and  then  performs 
a linear  interpolation  to  obtain  the  coordinates  for  the  intersection. 

When  the  user  opens  the  VRML  display,  the  3D-surface  and  a control  panel  appear.  Figure  4 shows 
the  look  of  the  VRML  display  when  the  VRML  browser  CosmoPlayer  is  used.  CosmoPlayer  can  be 
downloaded  from  http://cosmosoftware.com.  We  defined  another  VRML  prototype,  PROTO  Dash  to 
deal  with  the  operation  of  the  control  panel.  On  the  control  panel,  there  are  two  lines  of  buttons.  When 
the  user  clicks  the  X or  Y button  on  the  top  line,  the  specified  cutting  plane  will  be  shown.  For  example, 
clicking  the  X button  causes  a bounding  box  to  appear  around  the  figure  and  produces  a cutting  plane 
that  will  move  perpendicular  to  the  X axis.  The  buttons  on  the  second  row  control  the  animation  similar 
to  those  on  a VCR:  rewind,  reverse,  stop,  play,  and  fast  forward.  When  the  user  clicks  the  VCR  buttons, 
the  cutting  plane  moves  in  sync  with  the  projected  intersection  curve,  which  is  displayed  on  opposite 
faces  of  the  bounding  box.  The  text  on  the  bottom  of  the  browser  shows  the  value  of  the  cutting  plane 
location. 
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Figure  4:  VRML  display  on  CosmoPlayer. 


Figure  5 shows  the  imaginary  part  of  Airy  function  |Ai'(z)|  with  a cutting  plane  in  the  Y direction. 
Figure  6 presents  the  real  part  of  Airy  function  Bi(z)  with  a cutting  plane  in  the  X direction.  These  and 
other  figures  may  be  found  at  URL  http://zing.ncsl.nist.gOv/~  wang/mathlib/newmi/  . 


Future  research 

The  initial  attempt  to  develop  interactive  VRML  visualizations  for  the  DLMF  has  been  very  successful. 
Other  DLMF  project  members  are  pleased  about  the  results  to  date  and  have  agreed  to  work  closely  with 
us  in  determining  how  and  where  visualization  can  best  be  used  in  the  DLMF.  We  are  already  developing 
more  efficient  contour  plotting  code  to  facilitate  the  expansion  of  the  cutting  planes  to  the  Z direction. 
Also,  we  are  experimenting  with  the  use  of  a slider  bar  to  give  the  user  more  control  over  the  location  of 
the  cutting  plane.  Another  area  where  more  research  is  needed  is  in  the  development  of  efficient  clipping 
routines.  This  will  especially  be  a problem  when  we  begin  looking  at  visualizations  of  functions  like  the 
Gamma  function  where  the  domain  of  computation  will  contain  holes.  Each  chapter  will  produce  new 
challenges,  but  the  hope  is  that  much  of  what  we  are  learning  now  can  be  easily  applied  to  creating 
suitable  visualizations  for  the  other  chapters. 
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Figure  5:  VRML  display  with  Y direction  cutting  plane. 


Figure  6:  VRML  display  with  X direction  cutting  plane. 
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